Conformal fixed point of SU(3) gauge theory with 12 fundamental fermions 
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We study the infrared properties of SU(3) gauge theory coupled to 12 massless Dirac fermions 
in the fundamental representation. The renormalized running coupling constant is calculated in 
the Twisted Polyakov loop scheme on the lattice. From the step-scaling analysis, we find that the 
infrared behavior of the theory is governed by a non-trivial fixed point. 
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In recent years, there have been growing interests in 
the phase structure of asymptotically-free gauge theories 
with various numbers of fermions (Nf) and their repre- 
sentations. In particular, it is important to gain knowl- 
edge of the critical Nf (denoted as Nf) where different 
infrared (IR) properties are distinguished. Below Nf 1 , 
the theory is confining and chiral symmetry is broken, 
while above it there exists a non-trivial infrared fixed 
point (IRFP) which renders the theory conformal at low 
energy. The range Nf < N f < Nf is called the con- 
formal window, where Nf is the maximal Nf for the 
presence of asymptotic freedom. 

In addition to its intrinsic field-theoretic interests, 
studies of (near) conformal nature of strongly coupled 
gauge theories are of great phcnomcnological significance. 
Such a theory is a candidate for the dynamical origin of 
the electroweak symmetry breaking. This has led to stud- 
ies of conformal dynamics for various theories, using ap- 
proaches such as the Schwinger-Dyson equation, the ex- 
act renormalization group, etc.. Among these methods, 
lattice gauge theory has been widely used as a reliable 
tool to perform first-principle computation of scheme- 
independent, universal quantities. 

In this article we report results of a lattice study of 
SU(3) gauge theory coupled to 12 massless fermions in 
the fundamental representation [H-@|. There have been 
some recent lattice investigations in this theory, and the 
results have been controversial hitherto. In Ref. Q, 
the running coupling constant was computed in the 
Schrodinger functional (SF) scheme 0], and, within error, 
exhibited scale-independent behavior. Other works 0, 0| 
also showed the evidence of the conformal behavior in this 
theory. On the other hand, studies of the mass scaling 
behavior [B| and the spectrum of the Dirac operator [f| 
showed contradictive evidences. 

The renormalized coupling constant itself is a scheme- 
dependent quantity. However, by identifying the rela- 



tions among different renormalization schemes as coor- 
dinate transformations in the theory space, it can be 
demonstrated that the existence of a fixed point is scheme 
independent, except the case that the transformation is 
singular ||. In view of this, the best method for the 
search of the IRFP is to look for zeros of the beta func- 
tion. This leads us to focus on the ratio between the 
step-scaling function Q and the input renormalized cou- 
pling constant. This ratio approaches one when the beta 
function approaches zero. To confirm the existence of an 
IRFP in an asymptotically-free gauge theory, it is impor- 
tant to demonstrate that the ratio is one at both ultravi- 
olet (UV) and IR regimes, while being obviously different 
from this value in between. 

We calculate the renormalized coupling constant in the 
Twisted Polyakov Loop (TPL) scheme [l0( with the stan- 
dard staggered dynamical fermions. This scheme has no 
0(a) discretization error, which is a great advantage in 
taking the continuum limit. Another advantage is the 
absence of zero modes thanks to the twisted boundary 
condition, which enables us to perform simulation with 
massless fermions. 

The TPL scheme was proposed in Ref. for SU(2) 
gauge theory Here we generalize it to the SU(3) case. 
To define the TPL scheme, we introduce twisted bound- 
ary condition for the link variables, U fi (x) 1 in x and y 
directions on the lattice, 



U^x + PL/a) = Q„Un(x)Cll, {u = x,y) 



(1) 



where L is the spatial extent and a is the lattice spac- 
ing. The twist matrix f2„ is chosen to satisfy r^fjt = I, 
(f^) 3 = I, Tr[fi M ] = 0, and %9, v = e l27r / 3 ft„f^ 
for a given [i and v(^ /i). In the present study, we 
choose to work with the twist matrix as in Ref. [ill ]. 
The gauge transformation of is defined as U^(r) —> 
A(r)E/ At (r)AT(r + A)j with A(r) satisfying the relation 
h(r+vL / a) = fi„A(r)fiJ for the consistency with Eq.(p}. 
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By requiring the gauge and translational invariance, 
the Polyakov loops in a twisted direction can be defined 

as 

P x (y, z, t) = Tr(fY[ U x (x = j, y, z, t)]n x e^ 3L ) . (2) 

3 

Then the renormalized coupling in the TPL scheme is 
defined by taking the ratio of Polyakov loop correlators 
in the twisted (P x ) and the untwisted (P z ) directions: 

i (Ey. z My,z,L/2a)p x (o,o,oy) 

5TPL *<E SlV i'.(^»^/2o)P J ,(0,0,0)t>' U 

The leading discretization error is 0(a 2 ) in this scheme. 
At tree level, this ratio of the Polyakov loops is propor- 
tional to the square of the bare coupling. The propor- 
tionality factor k is obtained by analytically calculating 
the one-gluon-exchange diagram. In the case of SU(3) 
gauge group, the value of k with Wilson plaquette gauge 
action is 

k latt -(a/L) ~ 0.03184 + 0.00453(a/L) 2 + 0(a 4 ). (4) 

The discretization effects in k are of 0(a 2 ) instead of 
0(a), as expected. 

To include fcrmions in the fundamental representa- 
tion with twisted boundary condition, we introduce the 
"smell" degrees of freedom [Tgj] to avoid inconsistency 
with translational invariance. The smell quantum num- 
ber is a copy of color, and its non-trivial effects only ap- 
pear at the boundary. It is incorporated by identifying 
the fermion field as a N c x N s matrix, ip^(x), where a 
is the color index and a is the smell index. The twisted 
boundary condition for the fermion field is then imposed 
as 

r a (x + PL/a) = e™/ 3 Sl a u %(n u )l a , (5) 

for v ~ x,y directions. The smell degrees of freedom 
introduced here can be considered as extra flavors. This 
means that the number of flavors we can study on the 
lattice is a multiple of N s (= N c = 3). Furthermore, since 
we use staggered fermions in our simulation, we have 
four tastes for each flavor. This enables us to perform 
the computation with Nf = 3 x 4 = 12 in SU(3) gauge 
theory with twisted boundary condition. 

First, we discuss the vacuum structure and the exis- 
tence of the true vacua in an SU(3) gauge theory in- 
volving masslcss fcrmions in the deconfining phase. The 
free energy of the pure-gauge sector contains 3 4 -fold de- 
generate classical minima at V ' ^ = exp(27ri0 A1 /3)I, where 
6/j = 0,1,2 for each space-time direction. We investi- 
gated the semi-classical free energy in SU(3) gauge the- 
ory with Nf = 12 up to one loop level, and found that the 
vacuum energy is independent of #1,2, and the vacua with 
both #3,4 being 1 or 2 have the lowest free energy, indicat- 
ing that those are the true vacua. For these "non-trivial 



vacua" , all the classical link variables in z and t directions 
contain non-trivial phases ~ exp(±27ri/3L), giving 
rise to factors exp(±27ri/3) in the Polyakov loops. 

In this work, we generate gauge configurations around 
the true vacua where the vacuum expectation values of 
the Polyakov loops in untwisted directions have non- 
trivial phases. We observe that the vacuum stays having 
the nontrivial phases after thermalization and no transi- 
tion to trivial vacuum occurs. On the other hand, there 
is transition between nontrivial vacua in the low /? region. 
We monitored the value of the renormalized coupling dur- 
ing the transition and found that the contribution is neg- 
ligible. 

Our numerical simulation and analysis are conducted 
by the following guiding principles: 

1. We generate data in a broad range of /3, with in- 
tervals that ,g^p L grows almost constantly in each 
interval. Thus the interval of j3 is large in high f3 
region while small in the low f3 region. Each data 
have similar accuracy (2 — 3%). 

2. We employ fit functions for /3 interpolation which 
reproduce the tree level result <7xpl — So 011 eacn 
lattice size in extremely high (3 region. 

3. We include only the data which ensure that sys- 
tematic errors in all the interpolations and extrap- 
olations are under control. 

These guiding principles ensures the stability of our fit 
results under the change of fit functions. This is because 
the fit functions do not favor any special region of the 
data when we interpolate our data in /3 or extrapolate 
to in (a/L) 2 from the point 1, and the effect of statisti- 
cal fluctuation in high f3 region where the growth of the 
coupling is small is reduced from the point 2, and the 
data L/a = 4 data, having the largest discretization ef- 
fect, which was used in our previous report fl3| is now 
dropped. 

Our numerical simulation is performed in the following 
setup. The gauge configurations are generated by the Hy- 
brid Monte Carlo algorithm, and we use the Wilson gauge 
and the staggered fermion actions. The simulations are 
carried out with lattice sizes L/a = 6, 8, 10, 12, 16 and 20 
at around twenty f3 values (/? = 6/g 2 where go is the bare 
coupling) in the range of 4.0 < (3 < 100. To reduce sta- 
tistical fluctuations, we generate 8,000-897,000 trajecto- 
ries for each L/a) combination, measure the Polyakov 
loops at every trajectory and bin the data by taking the 
autocorrelation into account. The integrated autocorre- 
lation time is about 400 trajectories on the largest lat- 
tices, while as small as 3-5 trajectories on the smallest 
lattices. Using the Jackknife method, typical statistical 
errors of correlators are 2 — 3%. The bootstrap analysis 
produces consistent results. 

In Fig. [TJ we show simulation results for the renormal- 
ized coupling as a function of 1//3 for each L/a. For the 
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FIG. 1: TPL coupling for each /3 and L/a 



purpose of /3-intcrpolation, we use the following form of 
fitting function: 

<&r h <fi,a/L) =6/P + Y% =1 C i (a/L)/l3 i +\ (6) 

where N is the degree of the polynomial and Ci(a/L) are 
the fitting parameters. In order to obtain the best fit 
functions in which the x 2 /d.o.f. ~ 1, we take N = 3 — 5 
depending on the lattice size. 

To investigate the evolution of the renormalized run- 
ning coupling, we employ step-scaling function with two 
lattice sizes scaled by a step-scaling parameter s. For 
each L/a in the set of smaller lattices, we find the 
value of /3 which produces a given value of the renor- 
malized coupling, u=g 2 -, pL (/3, a/L). Then, we measure 
the renormalized gauge coupling at that /3 on the larger 
lattice, Z(u,a/L;s)=g^p L (f3,a/sL)\ g 2 n{ p a/L)=u . The 
step-scaling function, a(s, u), is obtained by taking the 
continuum limit of a/L; s): 



a(s, u) 



lim S(u, a/L; s) 

a— s-0 



(7) 



In this study, we take s = 1.5, and denote cr(u) = 
er(s=1.5,u) in the rest of this letter for simplicity. The 
set of smaller lattice is taken to be L/a = 6, 8, 10, 12, and 
therefore, we need values of <?tpl f° r L/a = 9, 12, 15, 18 
to take the continuum limit in Eq. ([7|). For L/a = 9, 15 
and 18, we estimate values of c/tpl f° r a given ft by 
the linear interpolation in (a/L) 2 with using the data 
on the lattices L/a = {8,10}, {12,16} and {16,20}, re- 
spectively. To estimate the systematic error of this in- 
terpolations, we also performed the linear interpolation 
in a/L, and found that the difference between a/L and 
(a/L) 2 interpolations is negligible. 

In Fig. [2j we show an example of the continuum ex- 
trapolation for obtaining a(u) in the strong coupling 
region (u = 2.054). The procedure we use to derive 
the central value of cr(u) is the linear extrapolation in 
(a/L) 2 with four points; L/a = 6, 8, 10, 12 -> 9, 12, 15, 18. 
Note that, in this example, though each lattice data 
a/ L; s = 1.5) is larger than u, in the continuum 
limit, <r(u) is consistent with u. This indicates that it is 
very important to take the continuum limit carefully 
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FIG. 2: Continuum extrapolation for the case of input cou- 
pling u = 2.054. Several kinds of extrapolation using all or 
various subsets of the data points are plotted. 
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FIG. 3: The growth ratio a(u)/u as a function of u with statis- 
tical (solid) and total (dashed) errors. Two-loop perturbative 
value (dashed line) is also plotted for comparison. 



We perform the step-scaling procedure explained above 
in a wide range of u. As mentioned in the beginning of 
this article, the growth rate a(u)/u is a suitable quan- 
tity for the search of the IRFP. This growth rate be- 
comes one when there is a zero in the beta function. 
Figure [3] shows a(u)/u as a function of u with statis- 
tical (solid) and total (dashed) errors. The total error 
includes both the statistical and systematic errors. The 
statistical error is estimated by Jackknife method. We 
will explain our estimation of the systematic error in de- 
tail later. In the weak coupling regime, the result is con- 
sistent with perturbation theory. At u = 2.05, the central 
value of a(u)/u touches 1, demonstrating the existence of 
an IRFP. Though several data points are also shown for 
u > 2.05 in Fig. [31 we note that the continuum extrap- 
olation cannot be performed reliably since the theory is 
not asymptotically free in this region. 

We also study the running behavior of the renormal- 
ized coupling obtained from each 100 Jackknife ensem- 
bles. They all converge in the IR energy region, confirm- 
ing the existence of the IRFP at 



2.05 ±0.47 (stat. 



)t° f 3 (syst.) 



(8) 



Here, the Jackknife error of the running coupling at IR 
is used as a statistical error. We will discuss our esti- 
mation of the systematic error in more detail later. The 
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corresponding /3 value for each lattice size at u* = 2.05 
can be calculated from the j3 interpolation function in 
Eq. ©: ( / 9,L/o)=(6.38,6), (6.65,8), (6.82,10), (6.96,12), 
(6.98,16), (7.02,18), (7.05,20). These are the parameter 
sets which reproduce conformal physics. 

Since the value of the coupling constant at IRFP 
is scheme dependent, to compare with other schemes, 
we calculate the critical exponent of the beta function 
as a scheme-independent quantity. In the vicinity of 
the IRFP, the beta function can be approximated by 
/3(u) ~ — "f(u* — u), where 7 is the critical exponent. 
We estimate 7 from the slope of a(u)/u against u, and 
obtain s"T = 0.95 ± 0.03. This leads to 7 = 0.14+°;°?, 
where the error is only statistical. Furthermore, if we use 
the different linear continuum extrapolation for L/a = 
6, 8, 10 -> 9, 12, 15, we obtained 7 = 0.40 ± 0.17. This 
means our 7 has a large systematic error. Still, the value 
we obtained here is consistent with both the perturbativc 

results, 72-ioop ~ 0.36 and 7£^ loop <~ 0.28 as estimated 
using 2-loop and 4- loop (MS scheme), and the nonper- 
turbative result in the SF scheme, 7sf = 0.13 ± 0.03 Q, 
within the systematic error. We remark that the value of 
7 is quite sensitive to the variation of the slope. Precise 
determination of the critical exponent is an important fu- 
ture problem to identify the interacting conformal fixed 
point. Another scheme-independent quantity which is 
interesting to observe is the mass anomalous dimension 
[H, EH near the IRFP. We will report it in a forthcoming 
paper. 

Here, we briefly mention the numerical stability of 
our results. For this purpose, we perform another step- 
scaling analysis based on s = 2 with L/a = 6,8,10 — > 
12, 16, 20. The continuum limit is taken by linearly ex- 
trapolating these three points in (a/L) 2 . We find that 
the running behavior in s = 2 step scaling is consistent 
with that in the case of s = 1.5, and the IRFP is found 
at u* = 2.26 ± 0.64 (stat.). This is consolidating the 
existence of the IRFP in this theory. We also derived 
the critical exponent of the beta function, and obtained 
7 = 0.11 ± 0.11. This is also consistent with our main 
results with s = 1.5 shown in the previous paragraph. 

Finally, we explain the systematic error in our analy- 
sis. There are two dominant sources of the systematic 
error. One is from the ambiguity of choosing the fit 
ansatz for the f3- interpolation (Eq. (|6])). The system- 
atic error regarding this is estimated as follows: Once we 
find, in each L/a, an optimal number of N in Eq. ([5]) 
which gives a reasonable fit (which will be used for the 
derivation of central values), we change the number TV 
to N + 1, and repeat the same fitting procedure. The 
difference between two fit results is adopted as the sys- 
tematic ambiguity in each L/a. (Wc should mention 
here that the systematic ambiguity of /3-intcrpolation for 
L/a = 10 is not included in our estimation of the system- 
atic error since only N — 4 gives a reasonable fit.) The 



other dominant systematic error comes from the contin- 
uum extrapolation. In Fig. [2] we show the comparison 
of several types of continuum extrapolation. As the cen- 
tral value, we take the linear extrapolation in (a/L) 2 for 
L/a = 6,8,10,12 9,12,15,18, and to estimate the 
systematic error in this procedure we estimate the differ- 
ences between the central value and the linear extrapo- 
lation without the finest lattice L/a = 12. We estimate 
the total systematic error by adding these differences in 
two origins (five f3 interpolations and the continuum ex- 
trapolation) of systematic error in quadrature. Further- 
more we compare the central value with the other lin- 
ear continuum extrapolation without the coarsest lattice 
L/a = 6 and quadratic extrapolation with all the data 
at four L/a. All the values in the continuum limit agree 
within 1-a statistical errors. For both of two analyses, 
the running region of u becomes narrower and the sta- 
tistical error of u* becomes much larger. Figure [3] shows 
the TPL renormalized coupling has a small systematic 
error in the strong coupling region, and we conclude that 
the existence of the IRFP is stable in this analysis. 

To summarize, we have found solid evidence of the ex- 
istence of an IRFP in SU(3) gauge theory with 12 mass- 
less Dirac fermions in the fundamental representation by 
using the TPL scheme. The procedure we presented here 
should be straightforwardly applied to studies of differ- 
ent gauge groups and fcrmion representations. It is also 
interesting to consider deformation of the theory at the 
UV scale with non-zero fermion mass or four-fermion in- 
teractions (see, e.g. [IBj]) in view of phenomenological 
applications. 
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